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I.  INTRODUCTION 


Many  projectiles  used  by  the  Army  are  slender  bodies  of  revolution 
which  are  launched  at  high  spin  rates.  Magnus  forces  and  moments  are 
generated  by  the  distorted  boundary  layer  which  results  from  a spinning 
body  at  angle  of  yaw.  Recent  Army  interest  in  achieving  increased 
range  and  greater  payload  capacity  in  artillery  projectiles  has  led  to 
designs  with  long,  slender  ogives,  increased  projectile  length,  and 
boattailed  afterbodies.  These  designs  have  resulted  in  decreased  drag 
with  a resulting  increase  in  range;  however,  the  aerodynamic  stability 
of  these  shapes  is  less  than  more  conventional  designs.  This  means 
that  these  new  shapes  are  more  susceptible  to  a Magnus  induced 
instability.  The  Magnus  force  is  small  (Figure  1),  typically  1/10  to 
1/100  of  the  normal  force;  however,  its  effect  is  important  because  the 
Magnus  moment  acts  to  undamp  the  projectile  throughout  its  flight. 

Thus,  it  is  desirable  to  minimize  the  Magnus  moment  in  order  for  the 
projectile  to  '’ly  at  a small  average  angle  of  attack  and  achieve  the 
greatest  range  capability. 

Magnus  has  been  modeled  theoretically  as  resulting  from  spin 
induced  distortion  of  the  boundary  layer.  This  effect  is  illustrated 
schematically  in  Figure  2 where  a cross-sectional  view  of  a body  of 
revolution  is  shown.  The  body  is  at  angle  of  attack  as  indicated  by 
the  cross  flow  velocity.  It;  the  view  where  there  is  no  surface  spin, 
the  profile  of  the  edge  of  the  boundary  layer  is  s>Tnmetric  with 
respect  to  the  plane  of  the  angle  of  attack.  In  the  view  where  the 
surface  is  spinning,  the  profile  of  the  boundary  layer  is  asymmetric 
with  respect  to  the  plane  of  the  angle  of  attack--thus,  the  inviscid 
pressure  distribution  (which  responds  to  the  aerodynamic  shape  composed 
of  the  model  + boundary-layer-displacement-surface)  is  asymmetric  and 
yields  a net  side  force. 

The  U.S.  Army  Ballistic  Research  Laboratory  has  placed  increased 
emphasis  on  research  into  the  Magnus  effect.  In  this  paper,  recent 
results  of  an  effort  to  develop  a method  for  computing  Magnus  effects 
for  use  in  projectile  design  will  be  discussed  and  results  of  related 
experimental  studies  will  be  presented. 


II.  THEORETICAL  APPROACH 


Background 


Since  the  Magnus  effect  is  a viscous  phenomena,  computation  of  the 
boundary  layei  development  is  the  foundation  for  computations  of  the 
Magnus  force.  The  boundary  layer  we  are  considering  is  fully  three 
dimensional  with  the  added  complication  of  the  interaction  of  surface 
spin  with  the  cross  flow  velocity.  The  inviscid  flow  also  requires 
special  attention  since,  in  order  to  compute  the  Magnus  force,  the 
inviscid  flow  computation  technique  must  be  able  to  compute  the  three- 
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dimensional  flow  over  a body  + boundary-layer-disj)lacement-surface  with 
no  plane  of  symmetry. 

B.  Boundary  Layer  Computations 

I'he  basic  equations  defining  the  three  dimensional  compressible, 
turbulent  boundary  layer  flow  over  an  axisymmetric  body  of  revolution 
described  by  the  relation  r = r(x)  are  listed  below.  The  coordinate 
system  is  shown  in  Figure  3. 

C'ontinuity 


^ (rpG)  ^ (rpv)  . (pw)  = 0 


(1) 


X- Moment urn 


- r-  3u  ' 3u  w 3u  w^  3r,  , i*  — ^ 

" 37  " 37  " 7 3?  - - = - 33r  " ^ ^ ^ 


e , 3 


(2) 


t- Moment urn 


3w 


3w  w 3w  uw  3r 


1 'I’e  3 


rw « n un  \.An  ux  ■>  l O r 'jVi  - ^ , y ^ 


3x 


3y  r 3())  r 3x 


r 3(J)  3y  3y 


Energy 


- I-  3h  " Dh  w 3h,  - ^^e  w ^^e  , .3u,  ^ f3w’ 


3u 


c^w  9 rU  9h 


- PU-V'  - - pv  W^  pv'h'l. 


f4) 


where  v = v + - — ■■  and  the  bar  indicates  a time  averaged  quantity. 

D 

In  order  to  obtain  closure  for  this  system  of  equations,  the 
following  models  of  the  turbulence  terms  have  been  introduced: 

lurbulcnt  Shear  Stress 


- I u 'v  ' 


pv  w = p 


'■3y-' 


3w  ^ 

1 


G 


r (-7^)' 
' 'Sy'' 


where  t is  introduced  as  the  turbulent  eddy  viscosity  and  the  mixing 
length,  I = .09  ! tanh  ( (.  1/ . 09)  (y/6) ] . Van  Driest  damping  is  used  to 
account  for  the  effect  of  the  laminar  sublayer. 
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Turbulent  Heat  Transfer 


pVTT  = 

c 8y 
P 


The  turbulent  Prcindtl  number  is  introduced  as 


Pr^  = c e/k^  = 0.90 
t P t 

The  numerical  technique  used  to  solve  these  equations  is  an 
implicit  technique  developed  by  Dwyer^’^  that  takes  into  consideration 
the  changes  in  direction  of  the  cross  flow  velocity  that  occurs  on  the 
side  of  the  model  where  the  inviscid  cross  flow  opposes  surface  spin. 
This  technique  correctly  models  the  cross  flow  convection  process 
occurring  within  the  boundary  layer.  In  order  to  improve  the  speed 
and  accuracy  of  the  numerical  solution,  several  coordinate  transforma- 
tions are  employed: 


Mangier  transformation  of  axisymmetric  growth, 


C = / r^  dx  ; 

o 

Blasius  type  transformation  of  normal  growth, 

"i  'I 

P P u y 

^00  OO  00  *•'  p 

J ^ rdy  ; and 
o 


P 2p  C 


coordinate  stretching  to  allow  closer  grid  spacing  near  the  wall^, 

n^  = 100  (1.5  exp  [ (j-1) (1/60) (1/.05) ] - 1)/(1.5  exp  (1/.05)  - 1) 


where  j = 1,  2,  3,  61 


1.  E.  A.  Duyer,  "Three  Dimensional  Flow  Studies  Over  a Spinning  Cone 
at  Angle  of  Attaakt"  BEL  Contract  Report  No.  127,  February  1974, 

U.S.  Army  Ballistia  Researah  Laboratories,  Aberdeen  Proving  Ground, 
Maryland.  AD  774796. 

2.  H.  A.  Dwyer  and  B.  R.  Sanders,  "Magnus  Foraes  on  Spinning  Supersonic 
Cones.  Part  I:  The  Boundary  Layer, " BRL  Contract  Report  No.  248, 

July  1975,  U.S.  Army  Ballistic  Research  Laboratories,  Aberdeen 
Proving  Ground,  Maryland.  AD  A012518.  Also  AIAA  Journal,  Vol.  14, 
No.  4,  April  1976,  p.  498. 

3.  F.  G.  Blottner,  "Variable  Grid  Scheme  Applied  to  Turbulent  Boundary 
Layers, " Journal  of  Computer  Methods  in  Applied  Mechanics  and 
Engineering,  1975. 
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In  computing  the  boundary  layer  development,  the  effect  of 
turbulence  is  turned  on  gradually  over  three  longitudinal  stejis.  The 
computation  grid  in  the  azimuthal  plane  is  in  10  degree  increments. 
Three  iterations  are  performed  at  each  station  for  turbulent  computa- 
tions. For  comparison  with  experiment,  the  location  of  boundary  layer 
transition  is  fixed  at  the  location  of  the  boundary  layer  trip  on  the 
experimental  model. 

C.  Three  lUmensional  Displacement  Surface 

The  three  dimensional  displacement  surface  is  not  merely  the 
vector  sum  of  the  longitudinal  and  circumferential  components  of  the 
boundary- layer  displacement  thicknesses.  Instead,  tlie  differential 
ecpiation  derived  by  Moore‘S; 


— P u r 
■X  e e 


Tl.  - ^ 


) 1 + f p w 
X 3((>  ' c c 


31) 


6^11  , 0 


(5) 
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mu.st  he  solved  for  ' , the  three  dimensional  boundary- 1 aver  displace- 

ment thickness  where 
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d>-  and  t 
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! n 

o 


j Vv 

c c 


dy 


Pw)*er  ha*^  shown  that  equation  (5J  is  of  the  general  form 


* ★ 


(bl 


★ -k 

I'he  i ncul.iri  t\‘  in  ' . and  at  x = 0 can  be  avoided  bv  starting  the 

■)!'  X 

computations  at  a small,  finite  value  of  x and  computing  approximate 
starting  conditions.  i.quation  Ibj  can  t lien  be  solveii  as  an  ordinar\- 
differential  ecpiation,  jiroviding  the  differencing  in  the  circumferen- 
tial direction  is  carried  out  from  : =0  to  180°  and  from  t = d fo  -180° 
in  order  to  obey  the  zone  of  influence  defined  by 


■/.  F.  N.  "DisfLcuement  Effect  of  a Three -?imenc,ional  iiour^Lir_- 

La’jey’,  " TN  E7Eo,  June  10b2. 

.5.  //.  A.  ['fir.fer,  "Methods  for  Corrqvitinij  Mafjnuo  Effcotn  on  Artillcru 

Product  ilea,"  BRL  doyitraet  Report  ,Vo.  S29,  Januaru  2977,  I'.F.  Arvyu 
Ballintic  Reneareh  Laboratory,  Abey\ieen  Py'ovituj  h'oimd,  Mai'nlarui. 
?h77,n. 


dx 

3*  = 


An  example  of  computed  values  of  6 jp  for  the  SOC  model  is  shown 

in  Figure  4.  The  effect  of  spin  is  to  decrease  the  thickness  of  the 
boundary  layer  where  surface  spin  and  inviscid  crossflow  are  in  the 
same  direction  and  to  increase  the  thickness  of  the  boundary  layer 
where  surface  spin  and  inviscid  crossflow  oppose. 

D.  Inviscid  Computations 

The  development  of  a numerical  technique  for  computing  the  three 
dimensional  inviscid  flow  field  over  a yawed,  pointed,  body  in  super- 
sonic flow  was  a very  important  step  in  the  development  of  a capability 
for  computing  Magnus  effects^.  The  program  uses  MacCormack ' s^  "shock 
capturing"  numerical  technique.  This  is  a second  order  accurate  scheme 
that  uses  a predictor-corrector  technique  to  solve  the  equations  of 
motion  in  an  implicit  marching  scheme.  The  unique  feature  of  the 
program  developed  by  Sanders  for  the  Magnus  problem  is  that  the  flow 
field  is  computed  about  an  axisymmetric  model  plus  displacement  surface 
which,  due  to  the  distortion  of  the  boundary  layer  by  surface  spin,  has 
no  plane  of  symmetry. 

E.  Sequence  of  Computations 

The  sequence  of  computations  which  must  be  run  in  order  to  compute 
Magnus  effects  is  indicated  in  Figure  5.  Each  block  indicates  a 
separate  computer  program  along  with  its  required  input  information  and 
the  output.  The  two  main  programs  are  outlined  in  asterisks. 

In  order  to  start  the  computations,  an  initial  plane  of  profile 
data  at  the  tip  of  the  model  and  inviscid  flow  boundary  conditions  are 
required.  The  initial  profile  data  are  computed  for  the  limiting 
condition  of  the  conical  tip  of  the  model.  The  inviscid  flow 
boundary  conditions  are  computed  using  the  program  developed  by 
Sanders®.  An  example  showing  the  computed  wall  pressure  distribution 
compared  to  experimental  data  for  several  azimuthal  stations  is  shown 
in  Figure  6.  Also  shown  is  the  outline  of  the  model  shape,  a six 

6.  B.  R.  Sanders,  "Three-Dimensional,  Steady,  Inviscid  Flow  Field 
Calculations  With  Application  to  the  Magnus  Problem, " PhD 
Dissertation,  University  of  California,  Davis,  California,  May 
1974. 

7.  R.  W.  MacCormaak,  "Numerical  Solution  of  the  Interaction  of  a Shock 
Wave  With  a Laminar  Boundary  Layer, " Proceedings  of  the  Inter- 
national Conference  on  Numerical  Methods  in  Fluid  Th^namics,  Lecture 
Notes  in  Physics,  Vot.  8,  Maurice  Holt,  ed. . Springer-Verlig , 1971. 
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caliber  long  secant-ogiye-cylimler.  As  seen  in  the  figure,  excellent 
agreement  is  obtained  between  theory  and  experiment  for  the  moderate 
angle  of  attack  considered  here.  It  is  seen  that  the  boundary  layer 
experiences  a significant  history  of  favorable  and  adverse  pressure 
gradients  in  both  the  longitudinal  and  circumferential  directions. 

In  order  to  start  the  boundary  layer  computatioi  for  tlie  spinning 
model,  initial  profile  data  are  generated  for  the  limiting  case  of  the 
laminar  boundary  layer  at  the  tip  of  a non-spinning  cone.  These  data, 
along  with  the  outer  boundary  condition  of  the  inviscid  flow,  enable 
the  marching  technique  to  begin  for  specific  conditions  of  Mach 
number,  angle  of  attack,  wall  temperature,  spin  rate,  and  free  stream 
properties.  The  output  of  the  boundary  layer  program  consists  of  wall 
shear  and  centrifugal  pressure  gradient  contributions  to  the  Magnus 
effect  and  the  longitudinal  and  circumferential  components  of  the 
boundary  layer  displacement  surface  as  functions  of  longitudinal  and 
circumferential  position  over  the  entire  surface  of  the  model. 

The  output  of  the  boundary  layer  program  is  input  to  tlie  program 
which  solves  for  the  three  dimensional  boundary  layer  displacement 

•k 

thickness,  6 Input  data  for  tliis  program  are  in  the  surface 

coordinate  system  used  for  the  boundary  layer  computations.  The  output 

of  this  program  is  transformed  into  a cylindrical  coordinate  system  in 

order  to  facilitate  computation  of  the  inviscid  flow.  The  output 

consists  of  the  surface  coordinates  of  the  model  plus  S*-,,  as  well  as 

‘ .•>11 

the  local  derivatives  of  the  surface  coordinate  plus  in  the  axial 
and  circumferential  directions. 

The  final  step  is  the  computation  of  the  inviscid  flow  over  the 
newly  defined  body  which  is  of  completely  arbitrary  configuration  with 
no  plane  of  s>Tnmetry.  The  starting  conditions  consist  of  the  inviscid 
flow  field  for  the  as>Tnptotic  cone  tip  of  the  original  model.  I’itch 
and  yaw  plane  force  and  moment  aerodynamic  coefficients  are  the  final 
outputs  obtained. 


Boundary  l.ayer  Components  of  Magnus 


Due  to  spin  induced  asNinmetry  in  the  com]nited  velocity  profiles, 
three  contributions  to  the  Magnus  effect  arc  generated  within  and  at 
the  surface  of  the  spinning  model  which  are  in  addition  to  the  boundar>’ 
layer  displacement  effect  sensed  by  the  outer  inviscid  flow.  These 
components  are:  (1)  longitudinal  velocity  wall  shear,  = a ( 'u/ 

(■2)  circumferential  vclocitv  wall  shear,  i = p f'w/">y'l  anil 

0 ' y=() 

y p 

centrifugal  pressure  gradient,  Ap  = J p — dy.  For  a non-spinning 

o 

model,  the  net  contribution  of  each  of  these  components  would  be  zero. 


However,  due  to  the  asymmetry  induced  by  surface  spin,  a small  contri- 
bution to  a side  force  is  obtained.  The  relative  magnitudes  of  these 
components  of  the  Magnus  force  are  shown  in  Figure  7. 


III.  EXPERIMENTS 


A.  General 


The  purpose  of  the  experimental  studies  is  to  provide  data  that 
will  be  useful  in  evaluating  and  providing  guidance  for  the  development 
of  the  theoretical  effort.  The  experimental  studies  consisted  of: 

(1)  boundary  layer  profile  measurements;  (2)  optical  studies; 

(3)  Preston  tube  skin  friction  measurements;  (4)  strain-gage  balance 
force  measurements;  and  (5)  wall  static  pressure  measurements. 

Experimental  results  were  obtained  in  the  BRL  Supersonic  Wind 
Tunnel  No.  1 which  is  a continuous  flow,  flexible  nozzle  tunnel  capable 
of  Mach  numbers  1.5  to  5.0;  the  test  section  size  is  38  cm  high  by 
33  cm  wide.  Data  were  obtained  at  Mach  3.0,  a Reynolds  number 

of  7.6  X 10®  and  at  both  spinning  and  non-spinning  conditions. 

The  model  tested  was  the  six  caliber  secant-ogive-cylinder  with 
geometry  as  shown  in  the  pressure  distribution  plot  of  Figure  6.  Tlie 
model  was  equipped  with  a boundary  layer  trip  located  0.7  caliber  from 
the  nose  to  insure  a consistent  turbulent  flow  for  all  tests. 

B.  Boundary  Layer  Measurements 

Measurements  of  the  total  head  pressure  through  the  boundary  layer 
were  made  with  a flattened  impact  pressure  probe  0.15  cm  wide  by  0.015 
cm  high.  Data  were  obtained  on  the  cylindrical  section  at  3.33,  4.44, 
and  5.56  calibers  from  the  nose  and  azimuthal ly  around  the  model  in  30 
degree  increments  and  at  yaw  angles  from  0 to  6.3  degrees.  At  each 
position,  surveys  were  made  at  both  zero  RPM  and  20,000  RPM;  the  spin 
rate  of  20,000  RPM  corresponds  to  a pd/V  of  0.19.  The  impact  pressure 
probe  was  brought  from  outside  the  boundary  layer  down  to,  and  contact- 
ing, the  model  surface  for  the  zero  RPM  case;  for  the  20,000  RPM  spin 
rate  the  probe  was  brought  down  to  within  approximately  0.01  cm  from 
the  surface.  The  probe  axis  was  aligned  longitudinally  with  the  model 
axis.  Some  uncertainty  is  inherent  in  the  profile  data  due  to  the 
probe  not  being  aligned  with  the  local  flow  direction  within  the 
boundary  layer.  This  uncertainty  would  be  greatest  near  the  surface  of 
the  model  and  at  longitudinal  stations  near  the  forward  portion  of  the 
model.  However,  the  large  gradients  present  in  a turbulent  boundary 
layer  would  confine  the  greatest  effect  of  flow  angularity  to  a very 
small  region  near  the  surface  which  cannot  be  probed  accurately  using  a 
total  head  probe.  Also,  these  measurements  were  obtained  at  small  angles 
of  attack.  Local  Mach  numbers  within  the  boundary  layer  were  deter- 


mined  from  the  Rayleigh  pitot  formula  assuming  a constant  static- 
pressure  through  the  boundary  layer.  The  boundary  layer  data  for  this 
paper  uses  measured  values  of  wall  static  pressure. 

Wall  shear  stress  was  obtained  for  the  non-spinning  model  using 
the  Preston  tube  technique.  The  Preston  tube  is  a circular,  total-head 
probe  mounted  flush  with  the  model  surface  and  sitcti  to  lie  within  the 
logarithmic  portion  of  the  law-of-the-wal 1 velucit)'  profile.  The  wall 
shear  stress  was  computed  from  the  measured  impact  and  surface  pressures 
and  using  the  correlation  relations  of  reference  K. 

V'apor  screen  flow  visualization  experiments  using  the  technique 
described  in  reference  9 were  performed.  The  pur]iose  of  this  jihase  was 
to  obtain  a better  understanding  of  the  surrounding  flow  field  and  in 
particular  determine  the  presence  of  vortices  both  imbedded  in  tlie 
boundary  layer  and  separated  from  the  model  surface,  \apor  screen 
pictures  were  obtained  at  each  caliber  position  for  angles  of  yaw 
of  2,  4,  5,  6,  and  10  degrees  both  with  and  without  spin. 

C.  Force  Measurements 

Measurements  of  Magnus  and  normal  force  were  obtained  using  the 
strain-gage  balance  technique.  The  model  is  free  to  rotate  on  intern- 
ally mounted  bearings  and  equipped  with  a single  row  of  turbine  blades 
so  that  spin-up  can  be  achieved.  Turbine  air  is  supplied  to  the  model 
through  the  sting  and  the  model  is  brought  up  to  speeds  as  high  as 
40,000  RPM.  The  turbine  air  is  then  cut  off  and  data  are  acquired  as 
the  model  spin  rate  decays. 


IV.  DISCUSSION  OF  TUT  TUF.ORFTICAL  AND  EXPFRIMIiNTAL  RESULTS 

A.  Boundary  Layer  Characteristics,  a < 4.2  degrees 

A comparison  of  theoretical  and  experimental  velocit)'  profiles  is 
shown  in  Figure  8 for  2 degrees  angle  of  attack  and  zero  spin  rate. 

The  agreement  is  considered  to  be  good  and  is  actually  comparable  to 
that  obtained  for  supersonic  two-dimensional  flow  measurements.  On  the 
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leeward  side,  near  <j>  = 180°,  the  theoretical  velocities  are  greater 
than  the  experimental  measurements  near  the  surface  and  the  theoretical 
velocities  are  smaller  than  the  experimental  velocities  near  the  edge 
of  the  boundary  layer.  On  the  windward  side,  near  i}>  = 0°,  the  situa- 
tion is  reversed.  These  differences  in  profile  shape  will  give  compen- 
sating effects  when  computing  integral  parameters.  Figure  9 compares 
theoretical  and  experimental  velocity  profiles  at  4.2  degrees  angle  of 
attack.  The  differences  between  theoretical  and  experimental  profile 
shape  are  similar  to  those  at  a = 2°  except  that  they  are  more  pro- 
nounced; also,  the  variation  in  boundary  layer  thickness  from  the  wind- 
ward to  the  leeward  side  is  greater  for  the  4.2°  case. 

The  effect  of  spin  on  velocity  profiles  is  shown  in  Figure  10 
where  profiles  on  the  side  where  cross  flow  and  spin  are  in  the  same 
direction  (<(>  = 0-180)  are  compared  with  profiles  on  the  opposite  side 
of  the  model  (ip  = 180-360)  where  the  cross  flow  is  in  the  opposite 
direction  of  spin.  On  the  windward  side  of  the  model  = 0-90),  there 
is  almost  no  measurable  effect  of  spin.  On  the  leeward  side  ($  = 120 
vs  240  and  150  vs  210),  the  profile  shapes  differ  substantially.  The 
effect  of  cross  flow  in  opposition  to  model  surface  rotation  (e.g. 

(j)  = 210)  is  to  decrease  the  fullness  of  the  profile  which,  of  course, 
will  result  in  a larger  displacement  thickness.  It  is  of  interest  to 
note  that  the  primary  effect  of  spin  is  to  change  the  profile  shape 
rather  than  to  change  the  total  thickness. 

Values  for  the  longitudinal  component  of  displacement  thickness 
are  compared  in  Figure  11.  The  agreement  between  theory  and  experiment 
is  generally  good;  however,  it  is  seen  at  the  forward  station  (Z/D  = 
3.33)  that  theoretical  thicknesses  are  slightly  greater  than  experiment 
(4>  0)  and  at  the  aft  station  theoretical  thicknesses  are  slightly 

smaller  than  experiment.  This  situation  indicates  that  the  boundary 
layer  actually  grows  at  a faster  rate  than  predicted  by  theory;  however, 
this  is  not  particularly  surprising  since  the  turbulence  model  did  not 
provide  for  any  adjustment  as  a function  of  pressure  gradient.  The 
effect  of  spin  on  displacement  thickness  6 can  be  seen  in  Figure  12 

where  the  increment  of  6 ^ due  to  spin  is  plotted  on  an  expanded  scale 

for  a = 4.2  degrees.  The  effect  on  displacement  thickness  is  seen  to 
be  significant  only  in  the  vicinity  of  the  leeward  side  ($  = 180°). 

The  agreement  between  theory  and  experiment  is  encouraging  evidence 
that  the  numerical  technique  accurately  models  the  effect  of  surface 
spin. 


Measured  values  for  skin  friction  coefficient  obtained  using  the 
Preston  tube  technique  are  compared  to  theory  in  Figure  13.  The  skin 
friction  coefficient  is  referenced  to  free  stream  static  properties 
upstream  of  the  model  rather  than  the  more  conventional  approach  of 
using  local  properties  at  the  edge  of  the  boundary  layer.  The  agree- 
ment indicated  is  within  ± lO^o.  This  is  considered  quite  good  since 
the  Preston  tube  is  expected  to  yield  an  accuracy  of  ± 10°4  for  two 
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Jimensional  flat  plate  iioundary  layer  flow.  The  use  of  tlie  Preston 
tube  to  obtain  measurements  in  a three  dimensional  boundary  layer  flow 
using  two  dimensional  calibration  data  must  be  regarded  as  speculative 
and  mainly  of  qualitative  interest. 

B . Boundary  Layer  Characteristics,  a > 4.2  degrees 

,\'o  boundary  layer  calculations  have  been  made  for  angles  of  attack 
greater  than  4.3  degrees.  [Experimental  data  at  larger  angles  of  attack 
shown  in  Figures  14,  1.3,  and  16  illustrate  difficulties  that  may  be 
encountered  in  theoretical  calculations.  A rejiresentati ve  sample  of 
the  vapor  screen  data  is  sliown  in  Figure  14  for  t = 10°.  Tlie  vapor 
screen  is  positioned  3.3  calibers  behind  the  nose  and  normal  to  the 
model  a.xis  of  s\'mmetry.  A svmmetrical  pair  of  separated  tiody  vortices 
are  visible  as  dark  regions  on  the  leeward  side.  Oata  obtained  at 
6.3  degrees  show  tlie  vortices  to  be  formed  but  they  arc  embedded  within 
the  boundary  layer.  Although  the  present  interest  is  at  lower  angles 
of  attack.  Figure  14  is  presented  here  for  clarity  purposes. 

Velocity  profiles  at  i = 6.34°  are  shown  in  Figure  15.  On  the 
windward  side  of  the  model,  the  profiles  seem  well  behaved;  however, 
on  the  leeward  side  U = 120-240),  jirofiles  differ  drasticall)’  with 
each  other.  For  example,  at  t = 180°  the  profile  shows  a fullness 
typical  of  profiles  seen  in  a favorable  pressure  gradient;  at  : = 210°, 
the  profile  is  less  full  and  seems  to  be  tending  toward  separation, 
typical  of  profiles  in  an  adverse  pressure  gradient.  Tlie  di  sji  1 acement 
thickness  at  = 210°  is  actually  greater  than  that  at  v'  = 1B0°  even 
though  the  total  thickness  is  smaller  than  that  at  >p  = 180°.  Another 
characteristic,  which  at  first  appears  abnormal,  is  the  large  differ- 
ence in  boundary  layer  thicknesses  at  <J>  = 210  and  240  degrees.  Dis- 
placement thicknesses  at  a = 6.34°  are  shown  in  Figure  16.  The  signi- 
ficant difference  from  the  4.2°  case  is  the  dip  or  decrease  in  ■ near 
: = 180°.  It  is  believed  that  the  above  phenomena  are  caused  by  the 
existence  of  longitudinal  separation  t\"pe  vortices  which  are  beginning 
to  develop.  Such  vortices  could  create  local  areas  of  favorable  and 
.idverse  pressure  gradients  that  would  cause  the  comple.x  flows 
illustrated  b>'  Figures  15  and  16. 


i’.  Magnus  Ch; 


One  criterion  which  can  be  used  to  gage  the  success  of  the  theorv 
is  the  accuracy  with  which  Magnus  forces  and  moments  can  he  jiredicted. 
Figure  17  is  a comparison  of  theoretical  and  experimental  Magnus 
forces;  in  addition,  the  magnitude  of  the  four  com]ionents  of  the  Magnus 
force  are  shown.  The  contributions  of  t,  and  Ap  oppose  and  are  of 

com|)ar;ible  magnitude,  while  the  contribution  of  t ^ is  minimal.  The 

arithmetic  sum  of  the  three  boundary  layer  components  is  indicated  b>- 
^YBL  ~ ^x  * computed  Magnus  force  showai  here  is  the 
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* 

arithmetic  sum  of  the  contributions  due  to  x , t , Ap,  and  6 It  is 

X ^ 3D 

worthwhile  to  emphasize  that  this  marks  the  first  time  that  computations 
of  the  Magnus  effect  have  been  carried  out  in  a conceptually  "exact" 
manner  for  the  turbulent  boundary  layer  on  a realistic  projectile  con- 
figuration. The  agreement  between  the  calculated  Magnus  force  and 
experimental  strain-gage  balance  force  measurements  is  extremely  good. 
The  theoretical  computations  accurately  reproduce  the  nonlinear  trend 
of  Magnus  with  angle  of  attack.  This  nonlinear  behavior  is  due  primari- 
ly to  the  increasing  dominance  of  the  contribution  of  ^ . as  the  angle 
of  attack  increases. 


V.  CONCLUDING  REMARKS 

A combined  theoretical-experimental  study  of  the  Magnus  effect  on 
yawed,  spinning  projectiles  has  been  discussed.  The  overall  objective 
of  this  effort  is  to  develop  a method  for  computing  Magnus  effects  that 
could  be  used  in  the  design  of  artillery  projectiles.  Numerical  tech- 
niques have  been  developed  for  computing:  (1)  the  three-dimensional 

turbulent  boundary  layer  development  over  a yawed,  spinning  body  of 
revolution;  (2)  the  three-dimensional  boundary- layer  displacement 
surface  for  an  arbitrary  body  of  revolution;  and  (5)  the  three-dimen- 
sional inviscid  flow  field  over  a yawed,  pointed  body  of  completely 
general  configuration  with  no  plane  of  symmetry.  The  computations  have 
been  compared  to  experimental  measurements  of  Magnus  force  and  turbulent 
boundary  layer  profile  characteristics.  The  general  impression 
obtained  in  comparing  the  computations  to  experimental  data  is  that  the 
numerical  techniques  are  working  quite  well  and  yielding  very  impressive 
agreement  with  experimental  data.  The  results  for  Magnus  force  are 
considered  extremely  encouraging.  The  comparisons  with  detailed  profile 
characteristics  do  reveal  minor  differences.  It  is  believed  that 
additional  accuracy  can  be  achieved  by  improved  or  more  complex  turbu- 
lence modeling  which  accounts  for  both  favorable  and  adverse  pressure 
gradients.  Other  refinements  in  the  boundary  layer  computation  such  as 
correction  for  transverse  curvature  and  inclusion  of  boundary  region 
effects  should  be  incorporated. 

The  ability  to  predict  the  point  of  transition  may  be  a factor  in 
the  accuracy  of  the  theoretical  model.  The  location  of  transition  was 
fixed  by  the  use  of  boundary  layer  trips  for  all  the  experiments 
described;  therefore,  the  theoretical  results  do  not  reflect  any  error 
which  might  be  caused  in  predicting  the  location  of  the  transition  line 
around  the  model. 

The  higher  angle  of  attack  experiments  suggest  that  problems  will 
be  encountered  in  extending  the  theory  to  angles  of  attack  greater  than 
5 degrees. 
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Figure  1.  Magnus  and  Normal  Forces  on  a Spinning  Projectile 
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Figure  2.  Schematic  Illustration  of  Spin  Induced 
Boundary  Layer  Distortion 
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Figure  3.  Coordinate  System 
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Figure  5.  Sequence  of  Numerical  Computations 
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Figure  13.  Preston  Tube  Skin  Friction  Results,  SOC  Model 
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TURBULENT  BOUNDARY  LAYER 
SECANT-OGIVE  CYLINDER  MODEL 
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Figure  17.  Magnus  Force  Versus  Angle  of  Attack, 
Theory  Compared  With  Experiment 
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LIST  OF  SYMBOLS 


skin  friction  coefficient 
Cp  specific  heat  at  constant  pressure 

normal  force  coefficient 
Cy  Magnus  (side)  force  coefficient 

D diameter  of  base  of  model 

h static  enthalpy 

turbulent  conductivity 

8.  mixing  length 

p pressure 

P spin  r?te,  radians  per  second 

Pr^  turbulent  Prandtl  number,  c c/k^ 

t ’ p ' t 

r local  radius  of  model 

Re^  Reynolds  number  based  on  model  length 

u,v,w  velocities  in  boundary  layer  coordinates 
V velocity  along  model  trajectory 

X surface  coordinate  in  longitudinal  direction 

y,Y  coordinate  perpendicular  to  local  surface 

z cylindrical  coordinate  along  model  axis 

a angle  of  attack 

e turbulent  eddy  viscosity 

6 boundary  layer  thickness 

6*  boundary  layer  displacement  thickness 
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LIST  OF  SYMBOLS  (Continued) 

Ap  centrifugal  pressure  gradient  contribution  to  side  force 

n transformed  y coordinate 

u molecular  viscosity 

C transformed  x coordinate 

p density 

longitudinal  velocity  wall  shear  contribution  to  side  force 
circumferential  velocity  wall  shear  contribution  to  side  force 

$ 

1)  coordinate  in  circumferential  (azimuthal)  direction 

Subscripts 

e edge  of  boundary  layer 

w model  wall  conditions 

X quantity  in  x direction 

free  stream  reference  condition 
<})  quantity  in  azimuthal  direction 

Superscripts 

fluctuating  quantity 
time  averaged  quantity 
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